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Abstract 



We study the lowest energy i? of a semirelativistic system of identical massless 
bosons with Hamiltonian 



N , N 




i=l j>i=l 



We prove 



< 



E 



< 



where A = 2.33810741 and B = ^ = 2.3447779. The average of these bounds 
determines E with an error less than 0.15% for all N >2. 
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One of the advantages of studying the semirelativistic "spinless-Salpeter" 
Hamiltonian [1,2] is that it captures some aspects of a fuU relativistic treatment 
and at the same time aUows us to express the many-body problem in a tractable 
form. The principal source of mathematical difficulty is the kinetic-energy operator 
\/ w? -\- , which is defined in momentum space as a multiplicative operator [3] , 
and becomes, via the Fourier transform, a non-local operator in configuration 
space. We have earlier found energy bounds [4] for systems of N bosons in the 
case m > 0. In the nonrelativistic limit m — > oo the kinetic energy has the 
Schrodinger asymptotic form a/ -|- ~ m + Since the Schrodinger 
many-body harmonic-oscillator problem is exactly soluble [5-7], we were able to 
derive energy bounds that are asymptotically exact as m — > oo. The bounds were 
weakest in the ultrarelativistic limit m ^ 0. It is the purpose of this paper to 
present accurate bounds for this limiting case m = 0. The Hamiltonian for the 
system we study is given by 

N N 

i=l j>i=l 

We shall prove that the lowest energy E of this system satisfies the inequalities 

A{^N^{N -ify < E < B {^N^{N -ify , (2) 
where the coefficients A and B are given by 

/81 \ ^ 

A = Ai(0) = 2.33810741, and B = i — ] = 2.3447779. (3) 

The energy E of the A'"-body system is therefore determined by the average of 
the bounds in Eq. (2) for all couplings 7 > 0, and all N > 2, with an error less 
than 0.15%. In order to establish the energy bounds we must consider two fun- 
damental symmetries: translational invariance, and boson permutation symmetry. 
The Hamiltonian H includes the kinetic energy of the centre of mass. Therefore 
we choose a set of relative coordinates so that the kinetic energy of the centre of 
mass can be eliminated. The most convenient relative coordinates for our purposes 
are Jacobi coordinates defined in terms of the (column) vector [r] of individual- 
particle coordinates by an orthogonal matrix R. Thus we write [p] = R[r]. Since 
R is orthogonal, the conjugate momenta [tt] are given in terms of the individual 
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momenta [p] by the expression [tt] = The first of the Jacobi coordinates, pi, 

is proportional to the centre-of-mass variable, so that the elements of the first row 
are all equal to 1/y/N. The other two coordinates which we shall need to refer to 
specially are p2 and pj^ which are given explicitly in terms of the 'other set' by 



ri-r2 1 N-l 

The expression of the boson permutation-symmetry constraint in terms of Jacobi 
coordinates can be a source of complication. But we do not need to face this difficulty 
here: we simply exploit the 'reducing power' of the necessary boson symmetry to 
relate the A^-body problem to a scaled two-body problem. Let us assume that \& 
is a normalized boson wave function of the — 1 relative coordinates {pi\f^2- 
By Lemma (1) established in Ref. [4], we know that an operator acting on \E' 
with leading term tti, may be replaced by zero, even when the term is inside the 
kinetic-energy square root. This will be important later. We shall also use another 
important relation [4, Eq. (2.5)], namely 

(*,p^^r) = (*,p2^), 2<z<iV. (5) 

An arbitrary boson wave function is not necessarily symmetric in the {pi}fL2' ^ut 
Eq. (5) is generally true, and is very useful. A special boson wave function which 
certainly is symmetric in the {pi}^2 Gaussian function which also 

has another unique [8,9] and useful property, namely it factors into single- variable 
Gaussians ij^g as follows: 

N 3 / 2 \ 

*fl(P2,P3,---,Piv) = X1^9(Pi)' V'sW = (-j exp f — —\, a>0. (6) 

That has the correct boson symmetry follows immediately from the following 
identity valid for Jacobi relative coordinates 

N N 
i=2 l=i<i 

In momentum space the Gaussian transforms to a Gaussian by the three- 
dimensional Fourier transform ^3 as follows: 

4>9 = Mi^g), Mk) = (^) ' exp (^-Q . (8) 
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The lower energy bound is found by the following argument. We suppose that 
is the exact ground-state wave function for the A^-body system corresponding 
to energy E and, using the necessary boson symmetry, we write 

E = H^) = (^^', |iV^+ (^^^Tln - r-2p| . 

By employing (4) and (5) in succession, and noting that the lemma allows us to 
remove the operator tti from the square root, we arrive at the relation 

E = {'if^H'if) = (^^,!^a^^J^^% + a-fp%^'^f^ , where a = N{N-l). 

Thus the A'"-body energy E is bounded below by the lowest energy £^ of the 
one-body Hamiltonian 

_ , 1 / — o 2 -^3 i 2 

rt = a'^ VP + 0!7T* a'^r + a^p . 

But ^^3(7^) is the Hamiltonian for the Schr5dinger problem of a single particle 
moving in a linear potential r. Thus we find 

E > = A{ja'^)^ , where A = Ai(0) fti 2.33810741 

is the first zero of Airy's function, and is also exactly the bottom of the spectrum 
of + r in three dimensions. This establishes the lower energy bound. 

The upper bound is found by means of the Gaussian 'trial' function ^'^ dis- 
cussed above. We have 



That is to say 



TT J 2a 

with respect to the variational parameter a > we obtain 

1 

1 / 81 \ 3 

E < = B{^a'^)s , where B=( — ] 2.3447779. 



2tt ) 

This result establishes the upper bound. It is perhaps tempting to try to improve 
the upper bound by the use of a more fiexible trial function. However, it is not 
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trivially easy to accomplish this, and to keep the calculation and result simple, 
since we must use a translation-invariant boson function. 

It is interesting that there is a relationship between the problem discussed 
in this paper and the corresponding Schrodinger problem with a linear potential 
and m > 0. To be more precise, if we consider the nonrelativistic problem with 
Hamiltonian H given by 

^ = E|^+ E ^i^^-^.i' ^>0' (9) 

then we have shown [10, Eq. (4.16)] that the lowest energy E oi H is bounded by 
the inequalities 

^(^^'(^^-If^)' < E < B{N\N-lf^\ (10) 

where the coefficients A and B are exactly as given in Eq. (3) above. The two 
problems are brought more nearly into 'coincidence' if we set 2m = 1. The re- 
maining differences, involving N and A , can then be understood if one notes that 
the (■^) factor must be associated with the 'potential terms' in each case, and, 
similarly, the A and 7 couplings must be made to 'correspond' by scaling, as the 
Fourier transformation, which relates the two systems, is applied. With the aid of 
such arguments, the 'almost equivalence' of the problems could eventually be used 
formally to extract our main result. However, we present these considerations here 
merely as a confirmation of our results; we prefer to develop more widely applicable 
direct approaches for the semirelativistic many-body problem itself, valid for all 
m > 0. 
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